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Q , Abstract: The constraints for evolution equations with some special forms of Lax pairs are first 

I investigated. We show by examples how the method is rooted in the classical literatures and how 

' the ignored constraints provide nontrivial solutions. Then we show, by the example of the KdV 

equation, how this special form of Lax pair may be found by the method of Wahlquist-Estabrook. 
At last we propose how to impose constraints for general Lax pairs including nonlinear ones. 
With the proposition the true Lax pairs and the fake ones can be distinguished easily. The 
linearity nature in integrable partial differential equations seems to have been revealed. 
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1 Introduction 



' In general, the last step of the inverse scattering method on the real line is to solve an 

1^ ■ integral equation such as Gelfand-Levitan-Marchenko equation. But only the reflection-less 

, potentials, which are solitons, make up of a closed system and therefore can be solved completely. 

I It had been observed that the reflection-less potential is some function of the eigenfunctions. 

With this observation the method of nonlinearization [2] of Lax pair was suggested. Then it 
was observed [8] that the constraint between the potential and the eigenfunctions may also 
I be regarded as symmetry constraint. Today the method of symmetry constraint has become a 

powerful tool for analyzing the solutions of integrable systems. By the method of nonlinearization 
of Lax pair or symmetry constraint, we will obtain a wider class of solutions much more than 
^ \ the solitons though the constraint maybe just be obtained by an observation on the soliton 

' and its corresponding eigenfunctions. In fact we will get algebraic-geometric solutions in most 



cases. The algebraic-geometric solutions are too complicated in practice, especially when we are 
only interested in the numerical integration of some initial-boundary problem for an integrable 
partial differential equations (PDE). Now it becomes more and more clear that a high-precision 
numerical integration of an integrable PDE should only cope with the constraint between the 
potential and the eigenfunctions. Yet we may expect those kinds of constraint will play a further 
role in both analytical and numerical applications. 

But the constraints are not so easy to find. The reason is that the form of the constraints 
may vary from one equation or hierarchy to another, and moreover for the same equation maybe 
there are several types of constraints, which may lead to completely different types of solutions. 
Now the method of symmetry constraint is still popular to get constraints, though a lot of useful 
constraints may be lost by it. In the first part of this paper we will analyze a wide class of 
Lax pairs, from which the constraints will arise naturally. The constraints provided here may 
be or not be the symmetry constraint. Therefore, sometimes the class of solutions may be 
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expanded. It must be point out that the special form of Lax pairs can also be gotten by the 
Wahlquist-Estabrook (WE) [TI] method. 

Sakovich [lOj had pointed out by example that Lax pairs with nontrivial spectral parameter 
may be fake ones. In fact the fake Lax pairs have long been puzzling. KaupfB] had advised a 
postulate to distinguish the fake Lax pairs from the true ones. But his postulate is not algorith- 
mic. And it also seems to be a little narrow in practice. Furthermore there are nonlinear Lax 
pairs though they may be linearized in most known examples. An algorithm for distinguishing 
the fake Lax pairs from the true ones should apply in nonlinear Lax pairs as well as the usual 
ones. In this paper we propose a functional equation, a solution to which will enable us to 
provide a nontrivial constraint for the integrable PDE. The constraint is simply a superposition 
of the solutions to the functional equation. Probably we may guess this superposition reflects 
more or less the linear nature of integrable PDEs. 

The paper is organized as follow. Section 2 propose a kind of natural constraints for a 
special form of Lax pairs. Several theorems there will be proved to guarantee the justice of 
the proposal. In addition several examples are provide in the section, where some interesting 
results may be found. Section 3 first shows by the example of the KdV equation how the WE 
method is capable of getting the special form of Lax pairs described in Section 2. Then a general 
requirement for constraining Lax pairs including nonlinear Lax pairs is proposed. At last the 
proposal is illustrated by two nonlinear Lax pairs of the KdV equation. Section 4 summarizes 
the main results. 

2 Constraints for evolution equations with a special form of Lax 
pair 

In the following paper, in order to avoid confusion by symbols, we define various symbols as 
fohows: ft = f = %; fx = f = /^"■^ = ^1 is the function obtained by applying 

operator P to function / = 1; Pseudo-operator d^^ is the inverse of operator d [Sj; g[u] denotes 
a differential polynomial of u. 

2.1 Observations about constraints for PDEs 

Numerical experiments have shown that most of the numerical integrations of initial bound- 
ary value problems (IBVPs) for PDEs are much more difficult than the IBVPs for ordinary 
differential equations(ODE). But for some PDEs the integrations may be done by only solv- 
ing a series of ODEs. This celebrated property should be considered as a kind of integrability 
for the PDE. As an example let us explain how the periodic KdV equation ut = 6uux + Uxxx 
is integrable in this sense. The stationary KdV equation is fn = 0, where fn = Ux and 
L = 5^ + 4u -|- 2uxd~^. It is well-known that /„ is a (2n-|-l)-order differential polynomial of u. 
So fn = may be written into a system of first-order ODEs for 2n-|- 1 variables u, u^^^ , • ' ' > u^"^""^ ■ 
Also we must know the time evolution for u, u^^^ , • • • , ti^^"^ . In fact 

4n(^) = -^(6nnW +n(=^)), Mod = 0, i = 0, 1, • • • , 2n. (1) 

at ax'- 

So we have 2 sets of ODEs for the KdV equation. One set governs the space evolution for 
u,u^^\--- ,u^'^"'\ And the other governs the time evolution for u,u^^\--- ,u^'^"'\ So for any 
(x, t) the values of u, u^^^ , • • • , u^^"^ may be gotten from the initial values uq, u^q ^ , • " " j ""o^"^ &t 
(xo,to) by a space evolution and a time evolution. In other words for some special kind of 
initial value problems the periodic KdV equation is solvable by only integrating some ODEs. It 
should be noted that the variables u,u^^\ ■ ■ ■ are only 'natural' variables, not 'canonical' 
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variables. The ODEs written for the 'canonical' variables are much simpler, see the following 
sections. 

Clearly from the point of view of numerical integration, the only crucial fact for the integra- 
tion of the periodic KdV equation is that /„ = is invariant under the KdV flow or in other 
words /„ = is an invariant manifold (IM) of the KdV equation. But for a given PDE its IMs 
are not so easy to give, especially when the required IMs must have some completeness. Luck- 
ily nonlinearization of Lax pair or symmetry constraint just provide a way to provide IMs for 
PDEs. Most popularly the symmetry constraint makes use of Lax pair in the inverse scattering 
transformation 

L^pi = Xitpi, (2) 
^V'i = P^i- (3) 
Still take the KdV equation as an example. The well-known nonlinearization of Lax pair 

LV'i = (a^ + = Ai^i, (4) 
= Pi^i = (40^ + 2,{ud + du))^i, (5) 



dt 



of the KdV equation is 



n 



co + '^Cii^h (6) 



i=l 

where tpi = tpi(x,t). With the constraint ([6]), Equation ^ and Equation ([5]) become two sets 
of ODEs. 

A less popular representation of Lax equation is 

h+i = Lh, (7) 

^/^ = Pfi. (8) 

for any i £ N. It is easy to verify that the compatibility condition for ([7]) and ([8]) is also the 
Lax equation in its operator form 

|i=|Ail. (9) 

So Equations ([7]) and ([8]) is equivalent to Equations ^ and ([5]). Sometimes a constraint for ([7]) 
and ([8|) is more convenient. 

There are several noticeable aspects for constraining the Lax equations. 

Our first observation is that if P is a differential operator then /n(x, 0) = is equivalent to 
fn{x,t) = 0. Let us take the KdV equation as an example to explain why this hold. The well- 
known P of the KdV equation is P = 45^ -|- 3{ud + du) . So /„ satisfies ft = ^fxxx + ^ufx + -iuxf- 
By analyzing the Taylor series of fn we know that if /„ is analytic and fn{x,0) = 0, then 
fn{x,t) = 0. Therefore, if P is a differential operator then ([7]) and ([8]) can be truncated by 
fn = 0. For the same reason the truncation may also be made by setting /„ = Y17=i ^ifi- 

Our second observation is that the number of equations in the truncated ([7]) and ([8]) is very 
close to but not enough for ([7]) and ([8]) to be 2 sets of ODEs. Also take the KdV equation as an 
example. To make it easy we take its natural Lax pair L = d"^ +'iu+2uxd~'^ , P = d{d'^+Qu). The 
number of equations in the truncated ([7]) and ([8]) including the KdV equation and the truncation 
condition is 2{n — 1) — = 2n — 1. But we need 2(n — 1) -|- 2 = 2n for /i, /2, • ■ ■ ; fn-i and 
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u in ([7]) and ([8]) to be ODE systems. So we need an additional relation between u and fi. This 
is just the nonlinearization of Lax pairs. For the example of the KdV equation the additional 
relation may be /i = Ux, which is just the usual way to get the stationary solutions. 

Our third observation is that for the periodic KdV equation the IMs mentioned above is 
dense. But for some other systems or another Lax pair, the relevant IMs may be not dense. For 
a nontrivial constraint of Lax equations it should be the minimum requirement that the number 
of the IMs increases as the truncation number n increases. 

The above three requirements rule out quite a lot of fake Lax pairs. But some true Lax 
pairs may also be ruled out. For example the famous short-pulse (SP) equation [9j Uxt = 
u + nn^ + l/2v?Uxx has a true Lax pair [T] 

L = {d-^ +Uxd-\x)d-\ (10) 

P = d-^ + -du^. (11) 
2 

Here P is not a differential operator. And furthermore SP equation is not an evolution equation. 
So the Lax pair (jlOp and (jlip does not fulfill our first requirement. Thereafter we will constrain 
ourselves in evolution equations. 

2.2 A special form of Lax pair and the relevant natural constraint 

We have mentioned that with the usual Lax pair ([5]) the KdV equation has constraint 
([6]). Except this famous constraint, we find that the KdV equation with the less popular Lax 
pair d?!) , dS]) also has another interesting constraint 

n 

u = cq+ ^ Ci+jfifj, (12) 

where fk = for k > n and Cfc = for k < n + 1. So for any n €z Z, fn = provides an IM for 
the KdV equation. 

At the first sight constraint ([6]) or (112p is strange. In this section we will investigate a special 
form of Lax pairs and give a natural constraint for them. Then constraints ([6]) and ()12p can be 
understood directly and easily. 

We start by the following simple lemma. 

Lemma 2.1 If P is differential operator and g[u] satisfies ^fl'['u] = Pg[u], then constraint 
di'^] — Y17=i ^i'^i '^'^ constraint g[u\ = ^11=1 ^^ifi > Y17=i ^ifi ~ *-* exisi, where ijji, fi are defined 

The proof is obvious. 

For a given Lax pair, it is not a easy task to find out a nontrivial constraint. However, there 
exist a natural constraint for Lax pairs in the following form (Form I): 

• L = + L^d~^, where is a differential operator and is a function; 

• P is a differential operator. 

Theorem 2.2 For a Lax pair in Form I, one possible constraint is = Y17=i ^i'^Pi; where 
Lipi = Xitpi. And the other constraint Y17=o'^ifi ~ ^> where fi = U is also possible. 
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Proof The first part. We should only prove -^[L^ = PL^ , because we already have ■§itpi = Pipi- 
The Lax equation is L = [P, L]. The negative part of the Lax equation is L^d~^ = {PL^)d~^ — 
L^d~^{P 1), where P is the operator conjugate of P. Then we get = —{Pl)d. So 

f F r> T F — f F r) T F — ^ 

— — = —{P 1) and ( — — )x = 0. So P 1 must be a constant. But we can redefine P by 
adding some constant to it such that P 1 = 0. Then we have — PL^ = 0. 

The second part. We will verify /j satisfies ([8]). /o satisfies ([8]), which has been proven 
in the first part. Suppose fi satisfies ([8]). Then ^fi+i = ^(Lfi) = {ft^)fi + L{^Ji) = 
{PL - LP)U + LPfi = PLU = Pfi+i. 

From the second part of the proof we know it is a useful method to nonlinearize the Lax 
system {(I7|), ([H])} by finding out a differential polynomial /o such that J^/o = Pfo- But this is 
not the only way, for example, see constraint (|12p . 



2.3 Some examples 

The following three examples are all Lax equations in Form L The first two examples, which 
are both of the KdV equation, show that different forms of Lax pair for the same equation may 
naturally lead to different kinds of solutions. The third example, which is the Ito equation, 
shows how a special form of Lax pair may be used by two ways to solve the same equation. 

Example The recursion operator of the KdV equation 

It is well-known that the recursion operator and the linearization operator form a natural 
Lax pair. Particularly to the KdV equation we have 

L = d'^ + iu + 2u^d-\ (13) 
P = d{d^ + 6u). (14) 

It is also fairly easy to verify that the Lax equation with the above L (I13p and P (jl4p is the 
KdV equation ut = Guu^ + u^xx- Obviously L (fTHj) and P ([HD is in Form I. Here = 
2ux- So a constraint Ux = Yl'^i^i exists, where ^pi satisfies (3^ + 4u + 2uxd~^)tpi = Xiipi and 
■§['4^1 = i'ixxx + (i{utpi)x- Another constraint is ^ Cj/j = 0, where /o = Ux, the famous finite gap 
constraint. 

This gives a simple explanation for why ([6]) holds. Note this classical constraint is the 
integration of symmetry constraint. Let (/>j = -^ipf, where ipi satisfies dH) and ([5]). It is easy to 
verify 

{d^ + 4n + 2uxd-^)^i = 4Xi(j)i, (15) 
= d{d^ + 6u)^i (16) 

So Ux = ai4)i = Y^^l^i <^i4^'^i i-e-, u = co + Ya=i "■i'^i is a proper constraint. 

How about (fT2]l ? Let fj = Y17=j "^ii^n+j-i- Then we can verify Lfj = /j+i and ^fj = Pfj, 
where 

L = d'^ + Au-2d-^Ux, (17) 
P = d^ + 6ud. (18) 

It is obvious ■§iiu — co) = P{u — cq). So u = CQ + ^aifi is a proper constraint, which is equivalent 
to (112]). 
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Example The soliton Lax pair of the KdV equation 

It is known that the KdV equation stiU has another Lax pair 

L = d + ud~^, (19) 
P = d^ + 3du, (20) 

which is also in Form L So n = Y27=i ^i'^Pi or n = ^27=1 ^ifi ^ proper constraint. 
We are extremely interested in the case I/"0 = 0. 

Theorem 2.3 (d + ud^^)"''^^0 = 0, n = 1,2, ■ ■ ■ , generates all n-soliton solitons of the KdV 
equation. 

Proof Recall the soliton solutions for the KdV equation [7] is 

n 

i=i 

First we will prove 

iV/' + uNi = -2kiN'i. 

Or 

LNl = -2kiNI + juu- 

Here the appearance of jnu is due to d~^, which is an indefinite integral operator. Then 

n 

L™n = 2^ [{-2Kjr + li{-2Kjr~' + ■■■ + 7m)c,iVj] . 

From the expression of L^u = we known that kj is completely determined by jj: 

{-2Kjr + 7ii-2Kjr'^ + . . . + ^„ = 0. (21) 

Altogether the n-solitons satisfy L"0 = 0. But n-soliton solutions have 2n free parameters and 
Z/"0 = can be written as 2n first-order ODEs. So the general solution of L^O = is the 
n-soliton solutions. 

Example Ito equation (Drinfeld-Sokolov II) 
The Ito equation [H [5] 

Ut = 3Va:, (22) 
Vt = {uv)x + Vxxx (23) 

has a Lax pair [4J 

L = d^ + ud + Ux + vd-^, (24) 
P = d{d'^ + u), (25) 
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which is in Form I. We will investigate the constraints L"0 = 0. The first few constraints of (j22p 
and (f23]l can be easily solved. 

The first constraint is LO = v = 0. By (|22p ut = 0. So u = uq{x). 

The second constraint is L^O = Vxxx + {uv)x + v f vdx = 0. By ([23l) vt = —v J vdx. Let us 
introduce a new variable (j) such that (f>x = v. Then = —4>x(t> = —■^{(l)'^)x- So = —^(j)'^+g{t). 
We can prove is independent on t. So 



= -\<^' + |. (26) 

The solution of ([26]) is = y^tanh(-^t+ci). cq is independent on x. But ci may be dependent 
on X. So the final result of ijj is 



We may prove 



^tanh(^t + ci(x)). (27) 



So 

2 2 

= (/'x. (30) 

is one solution of (j22p and (j23p . 

The third constraint is L^O = 0. By a tedious study of this constraint, we find if (j) satisfies 
the following equation 

(I)t = -l/2(t>^ + fl{t)(t> + f2, 

where 

and /2 is an arbitrary constant, then v and u expressed by 

V = 4>x, 

U = —{(t)t-(t)xxx-l)- (31) 
(Px 

is a solution. 

It seems very difficult to solve L"0 = 0, n > 3. So how to analyze the high-mode vibrations? 
First by Theorem 12.21 we can impose a constraint v = where Lipi = Xiipi and 

= Pipi- After a little modification of the form of the Lax equation, we get 
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Let ipi = Lp[. Then v = X^ILi ^i^'i- So we get 

<fi = ipi' + U(Pi, 

u = 2>v' . 

Multiplying (j32p by ai and summing over i we get 

n n ^ / n 



ai'Pi] 



(32) 
(33) 
(34) 



(35) 



Integrating 



we immediately get 



1=1 



1=1 



1=1 \i=l / 



+ 7- 



At first glance 7 is a function of t. But by (j34p we known 7(t) = 0. So 7 is a constant. Now u 
is solved as 



u 



^ n n 

-(^Oic^i)^ + v" + j - y^^aiXjipi 



i=l 



i=l 



(36) 



Totally Equation (j32p and Equation (j33p contain (n — l) + n = 2n — 1 differential equations. So 
it is not enough for them to be ODE systems. The usual way to overcome this problem is to 
introduce another constraint. But for the Ito equation there is no need to do such things. The 
reason is as follows. Obviously we have 



(37) 



By differentiate (j37p with respect to x once and twice we get another 2 sets of differential 
equations. Together with ([33|) and ([37|) we have a closed form for the t-evolution of ip'^, ip'l 
and (p'l' . This is enough for a numerical computation for the Ito equation. Note that Equation 
()32p may be regarded as 4n — 1 first-order differential constraints in the x-direction, i.e., there 
is still one arbitrary function ipiQ{x) among the 4n initial functions Lpi. 



2.4 Further generalizations 

Now we will generalize Form I a little. Anyway we have the following theorem: 

Theorem 2.4 If L = + '}2i=i fi'^^^dij where is a differential operator, and also P 
is a differential operator, then there a linear composition of fis = Yl^=i'^ifi such that 
= PL^ , i.e., a nonlinearization = Yl'iLi Pii^i exists, where "0 is eigenfunction Lip = Xi/j. 

We say L and P in Theorem 12.41 is of Form II. To prove Theorem 12.41 we need the following 
Lemma 12.51 and Corollary 12.61 

Lemma 2.5 Iffj:7=i f:{x)d-^ gi{x) = 0, then ^Hi f^{x)g^{y) = 0. 

Proof Tri=ifi{x)d-^9i{x) = is equivalent to TJl=^ f^{x)gf\x) = 0, A: = 0, 1, 2, • • • . So 

Er=l m9^{y) = Er=l f^{x)9^{x + {y- X)) = Er=l Er=0 y^{x)9f\x){y - Xf = 0. 
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Corollary 2.6 If T.7=i Mx)d-^ Bi{x) = YTi=iFi{x)d~^G^{x) and both {A,{x)] and { Gi{x)] 
are linearly independent, then Bi{x) is a linear composition of Gj{x) and Fi{x) is a linear 
composition of Aj{x). 

Theorem 12.41 guarantees that there is a natural constraint for a Lax pair in Form II. For 
Lax pairs in more comphcated form, such as L and P are matrix, we have not reach a general 
result like Theorem 12.21 But it seems always a good guess that a linear composition v of the 
components before satisfies i) = Pv. 

Example ZS-AKNS 

With the help of recursion operator the AKNS hierarchy can be expressed as a simple 
expression 

f _l f -d + 2qd-^r 2qd-^q \ 
i\ -2rd-'r d - 2rd'\ J ' 

^ -L-( -"^], (39) 



r J ^ \ ir 



where i is the imaginary unit. Here L is the recursion operator. So L and the linearization 
operator P form a natural Lax equation = [P,L]. The symmetry ~ (^^-^ satisfies 

<7i+i = Lai, 

—ai = Pai. 
dt 

But ( ^ ) is also a symmetry because if ( ) is a solution of (j38p and <\39\i then ' ^ 



kq 
r/k 



is also a solution, where /c G M. So we know 



d_ 

Jt \ -r 



Therefore one possible constraint is ( ^ | = ^^ajUj. Meanwhile we can also write the Lax 

^ ^ t=i 

equation as = A^' and = P^, where ^' is a 2 x 1 vector. So | J = fflj^'i is also 

\ ^ ^ i=i 

one possible constraint. 

For example AKNS n = 3 is the coupled KdV equation 



r J ^ V ^Q'^'^x - r, 



XXX 
XXX 



The corresponding linearization operator is 

Qrdq — Qqqx 

^ = v 

It is easy to verify 



^ i 6rrx 6qdr — ' ^^^^ 



= P 

-r 
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Then one possible constraint is 



-r 



i=l 



where ( J = ^'j satisfies L^i = Xi^i, = P^i and the expressions for L and P are ([38 
and (lioi). 



2.5 The equivalence of the two kinds of constraint 

We have introduced two kinds of constraints in Lemma 12.51 and thereafter. Now we will 
explain the relation between the two kinds of constraint. 
Recall the constraint of the first kind is 

n 

g[u]=Y,(^i^i- (41) 

i=l 

Constraint of the second kind is 

n 

9[u] = Yhfi. (42) 

i=l 

In general V'i and fi have not any relation. But with constraint (|4ip or (j42p they become linear 

n 

dependent in most cases. Suppose fi = Cijipj. Then from the relation /j+i = Lfi, we get 

j j 

So Cj+ij = XjCij. 

n 

Theorem 2.7 For a constraint of the first kind g[u] = ajipj satisfying that if i ^ j then 

i=l 

n n 

Aj 7^ Xj, there exist a constraint of the second kind g[u] = bjfj, fn+i = kjfj where bi and 

i=l i=l 

n 

ki is determined by Ci+ij = XjCij, aj = ^^ftjCy, A" = kjXj'^ with arbitrary cu ^ 0. 

i=l 

Proof Only to consider that if H ^li 7^ and Aj 7^ Xj then det(cjj) 7^ 0. 

n n 

Theorem 2.8 For a constraint of the second kind g[u] = '^^bifi, fn+i = ^^^i/i satisfying that 

i=l i=l 

n n 

A" = kjX^~^ has distinct root for X, there exist a constraint of the first kind g[u] = aitpi, 
j=i i=i 

n 

where Aj and Oj is determined by: Aj is the i-th root of X^ = ^^^A-'"^ , cu 7^ 0, Q+ij = XjCij, 

flj = bjCji, . 
j 
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3 Searching the special form of Lax pairs and constraining the 
nonUnear Lax pair systems 



It is well-known that the WE method is a powerful tool for searching the Lax pairs for 
PDEs. By the preceding examples we have known that one PDE can have completely different 
Lax pairs. So which Lax pair may be found by the WE method is heavily dependent on the 
original assumptions for the form of Lax pair. Furthermore the Lax pairs obtained by the 
WE method are often nonlinear ones, which must be linearized before further application. So 
very often the nonlinear Lax pairs are considered to be useless. In this section we will first 
demonstrate, by the example of the KdV equation, how to search the Lax pair in Form II. Then 
we will propose a functional equation for finding constraints for both linear or nonlinear Lax 
pair systems. 

Let us first demonstrate by the example of KdV equation how Lax pair (I19p . (I20p can be 
found by the usual WE method. The KdV equation is first written into 



Ux — Z, 
Zx = P, 

ut = Quz+px- (43) 
Substituting ()33|) to the zero-curvature equation 

Mt-Nx + [M, N] = 0, 

where M = M(u, z,p),N = N{u, z,p), we immediately get 

ztM, + ptMp + px{Mu - Np) + 6uzMu - zN^ - pN, + [M, N] = 0. 

Here zt, pt and px should be considered as independent variables. 

Remark: The KdV equation in the form of (j43p . in fact, has implied that M is only dependent 
on u. In other words, we are searching Lax pairs for the KdV equation whose M is only dependent 
on u. 

Then by the usual steps of the WE method we get 



M = Xi + uX2 + u'^X3, 



n2 



N = + (3n^ + p)X2 + (4u^ + 2up - z')X3 + z[Xi,X2] + u[Xi, [X^X^]] + y [^2, [^i, ^2]], 



where 



[^1,^3] = 0, 

[^2,^3] = 0, 

[Xi,X^] = 0, 

[Xu[Xi,[Xi,X2]]] + [X2,X^] = 0, 

[X2,[X2,[Xl,X2]]] = 0, 

^[Xi,[X2,[Xi,X2]]] + 3[Xi,X2] + [X2,[Xi,[Xi,X2]]] + [Xs,Xi] = 0. (44) 



We would first not consider (|44p as an open Lie structure, but rather a matrix equation. Then 
we will get linear Lax pair. The 2x2 realization of (j44p can be solved by Maple and there are 
12 solutions obtained by Maple. However, only one solution is worthy of notice. The solution is 

a^3;2,i = 0, 
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X2;l,2 = 0, 
X2;2,2 = X2;l,l, 
X3;l,2 = 0, 
X3;2,2 = X3-l^l, 
X2;2,l = -Xl-1,2~^ 



X4-2,l = Xi-2,1 {xi;l,l^ + a;i;2,2^ - 2xi;i,iXi;2,2 + 4xi;i,2Xi;2,l) , 



a;4;l,2 = Xi;i,2 + Xi;i,2Xi;2,2 - 2xi;i,iXi;i,2Xi;2,2 + 4xi;i,2 Xi-2,1, 

Xi-2,2 = X4-i^i - xi-i^i^ - 3 3;i;i,ixi;2,2^ + 3 xi;i,i^xi;2,2 - 4xi;i,ia;i;i,2Xi;2,i 

+xi;2,2^ + 4xi;i,2a;i;2,ixi;2,2, (45) 

where Xi-j^k is the (j, /c) element of matrix Xi. Thanks to ()45p . the equation 

*x = M^, (46) 

where = ( , / ' ( 1 is a 2-dimensional vector, is completely determined. The goal equation 

V H^{X, t) J 

is Lip = Xip. (/? may be the linear composition of tj: and </> 

V^ = /i(M)V' + /2(M)0. (47) 

Clearly ip satisfies a 2-order ODE. If we demand the 2-order ODE do not contain u^^\n > 1, 
then we obtain /2 = 0, fi = Const, x^-i^i = and X2;i,i = 0. Obviously, we may set fi = 1. 
Now (HSl) is reduced to 



Vxx - (Xl;l,l + Xi;2,2)v5x + iu- Xi;i,2Xi;2,l + Xi;i,iXi;2,2)V5 = 0. (48) 

Comparing Equation (I48p with the goal equation Lip = Xip, we immediately get 

^xx + Xnpx + {u + X2)ip = 0. (49) 



To write a Lax pair in Form II from ()49p . we only need to introduce the variable ^ = ipx- Now 
l9l) can be rewritten to 



Li^{d+{u + X2)d-^)i = -Xii. 

Now L here is in Form II. That P is also in Form II can be verified directly. 

With regard to the other 11 solutions of ()45p . four of them can not be reduced to 2-order 
ODE; seven of them can not be reduced to 2-order ODE whose coefficients are only functions 
of u. 

Now we turn to nonlinear Lax pairs. It has NEVER been demanded the Lax pair be linear. 
Then where is the linearity? We find the following functional equation is crucial in integrable 
systems. The functional equation is 

j^g[u,^,X]= A[u]g[u,i:,X], (50) 

where g[u, ip, X] denotes a function of finite variables including the original variables u^^^ , auxiliary 
variables ipi and complex parameters A. But A[u] is a linear differential operator which only 
involves functions of finite variables of n*^*^ . Of course, non-degenerate conditions -^g[u, ip,X] ^ 
and j^g[u, ^,X] ^ must be satisfied. The final constraint for the general Lax system is 

^<7KV',A]=0. (51) 
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It seems that the hnearity nature of Lax integrable systems is completely characterized by 
Equation ()50p and (I5ip . To see how these arguments works, let us still take the KdV equation 
as an example. 

Following WE we regard Equation (I44p as a Lie algebra constraint. Here we only cite the 
final equations for the pseudo-potentials: 

ox at 
^ = -2/8, — = z-2(n + 2A)y8, 

^ = X-u-yl ^ = 2zys - 2{u + 2A)(n + yf - X) - p. (52) 

Note the coefficients have been adjusted to fit the KdV equation here. Because (j52p is a nonlinear 
Lax pair, the subsequent step is usually to find some transformation to linearize (|52p . But 
Equation ([50j) and ([5T]) enable us to solve the KdV equation straightly. Our task is to find a 
function g and a linear operator A satisfying (jSOp . It is easy to check 



g = ae-'^y^ + (iu + -i, A = d'^ + Qud, a,/3,7GC, (53) 
is a solution of (f50]l . So 

n n 

+ ^.^ + = ^ + + ^ aie-2s'« = 0, (54) 

i=l i=l 



is a proper constraint, where y^a satisfies (|52|) with A = Aj. 

Given a nonlinear Lax pair, (j50p can always be proposed no matter whether the nonlinear 
Lax pair is linearizable. The difficulty is to give nontrivial g and A. 

Open Lie algebra (|44p has more than one solutions. Let us assume = 0. (It seems also 
to be not clear yet when X3 is a centre.) One 4-dimensional realization of (|44p is 



1 

, 2 q2v o„,2 



1 



[X2, X4] = /3X4 - ^X5, [X2, X5] = a'(3'Xi - 2a' X2 - fSX^, [X4, X5] = 2qXi + X2. 

If we set X4 = X4 — -^Xi, then X4 commute with any other vectors. If there is no centre then 
X4 = 0, i.e., X4 = j^Xi. Therefore, in fact we have a 3-dimensional realization of (I44p 

[X,,X2] = - f3Xi, [Xi,X,] = 2a^Xi + 0X2, [X2, X5] = a^/3^X^ - 2a^X2 - PX^. (55) 

Following WE we get the following nonlinear Lax pair 

dyi a 11. 7 _JiL 

— — = — 1 + au{2 ev^ e v^j, 

ox 7 a 

9yi 1 , . 2,0" I n 2 I ^ / , //N 4^ 

— — = h 4an + 2au (n + zau + + au jev^ 

Ot a 7 

(u + 2an^ — i/au' + an")e 

a 

dm ^ .c^ „^ 
5x 7' 

^ = ^ - + 2^2 + +u" ev^ - - + 2u2 - ^ + u" e v^. 56 
ai 7^ a va a 



n(— e — e ^), 
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By a definition of yi = + In a — In 7 and y2 = -^y2, we can simplify (j56p to 

dyi 1 ^ I— 

= j= + 2^ au — 2va^* cosn(?/i), 

ox a/ ol 

= ^( + + 2ti") - 2Va(- + + u") cosh(yi) - 2u sinh(yi), 

ot a 

— =usmh[yi), 

^ = JfL cosh(yi) + (- + 2v? + u") sinh(yi). (57) 



Still another substitution yi = e^^ , e = -y/a and y2 = 2y2 transforms (i57|l to a form more 
suitable for computation, 

dyi 1 2 

— = -en + (2en - -)yi - euy^, 

^ = _ 2eu^ + u'- eu" - e{\ - Au^ - 2u")m - (- + ^eu^ + u' + eu")yl 
ot e e 

dm ~ 1 
— — = uyi -u—, 
dx yi 

-TTT = + 2u H h n )yi - + 2n h n — . (58) 

Ot e e y\ 

Now the task is to give g and A. With the help of Maple we get the following solution 

, (yi - 1)^ ^2 
5 = Co + cin + C2 ; e 

i = ^3 _^ 6ua. (59) 



Then it becomes very clear that (j58ll and (j52]) must be equivalent. In fact, the transformation 
from (158)1 to (1521) is 



ys = TTT"^' 2/3 = ^ey2 + ln(yi - 1) - ^Int/i, A = (60) 
2e 1 — yi 2 2 4e^ 



4 Conclusions 

We have proposed a kind of natural constraints for evolution PDEs with a special kind of 
Lax pairs. By the method several examples have been studied in detail. At least two interesting 
results should be noticed. One is that KdV equation has a soliton-Lax pair, from which only 
soliton solutions can appear. The other is that the Ito equation can be constraint to a series of 
PDEs that are solvable by the method of characteristics. 

For a given Lax pair, it is always very difficult to solve the functional equation ()50p com- 
pletely. Even for the KdV equation only special solutions of (I50p have been obtained. It is 
necessary to point out that solutions (f53|) and (j59|) are both obtained by the classical separation 
of variables. More general ansatz about the solutions of (150p will lead to a set of too compli- 
cated equations to solve. This also explains why we seek the special kind of Lax pairs, of which 
the constraints are manifest. Equations (150p and (I5ip can act as a measure of integrability for 
evolution equations with Lax pairs. Any generalization of ([50P and ([5T|) seems so difficult. 
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